Anomalous lifetime distributions and topological traps in ordering dynamics 



X. Castell6*'\ R. Toivoncn*'!^^^ y. M. Egufluz\ J. Saramaki^, K. Kaski^ and M. San Miguel^-^ 
^Umdad de Fisica Interdisciphnar-IMEDEA (CSIC-UIB), E07122 Palma de Mallorca, Spain 
^Laboratory of Computational Engineering, Helsinki University of Technology, P.O. Box 9203, 02015 HUT, Finland* 

(Dated: February 1, 2008) 

We address the role of community structure of an interaction network in ordering dynamics, as 
well as associated forms of metastability. We consider the voter and AB model dynamics in a 
network model which mimics social interactions. The AB model includes an intermediate state 
between the two excluding options of the voter model. For the voter model we find dynamical 
metastable disordered states with a characteristic mean lifetime. However, for the AB dynamics we 
find a power law distribution of the lifetime of metastable states, so that the mean lifetime is not 
representative of the dynamics. These trapped metastable states, which can order at all time scales, 
originate in the mesoscopic network structure. 
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I. 



INTRODUCTION 



structural correlations. 



Statistical mechanics and complex network theory 
have been applied to different disciplines, ranging from 
biology to sociology. From this perspective, social sys- 
tems are modelled as a collection of agents, located at the 
nodes of a network, interacting through simple rules. So- 
cial networks of human interaction are structured into co- 
hesive groups [1], and increased knowledge of this struc- 
ture [2-4] has sparked the creation of new network models 
[3, 5-9]. These models allow us to study the effect of the 
structure of social interactions on the dynamics taking 
place on the networks, and on the associated collective 
phenomena emerging from the interactions among the 
agents. 

The mesoscopic structure of a social network, and in 
particular its community structure, has been found to 
influence dynamics taking place on it in ways that can- 
not be explained by global level statistics in several cases 
[4, 10, 11]. In this paper we address the role of such 
mesoscopic structure on ordering dynamics or consensus 
processes: the question is when the interaction of agents 
with several options leads to an ordered state with a sin- 
gle option (consensus) or when disordered states (pos- 
sibly metastable), with coexistent options prevail. We 
consider two dynamical models. The first one is the 
prototype voter model [12] whose dynamics in complex 
networks is known to be generally determined by global 
properties such as the effective network dimensionality 
[13]. Secondly, we consider the AB model [14] introduced 
to describe language competition, which gives a natural 
context for the community concept. These two dynami- 
cal models are studied in a class of networks [9] incorpo- 
rating nontrivial community structure which introduces 



II. TWO DYNAMICAL MODELS OF 
COMPETING OPTIONS 

The voter model [12] concerns the competition of two 
equivalent but excluding options A and B. The state of a 
node is updated by imitation of a randomly chosen neigh- 
bor. The AB model [14] includes a third non-excluding 
mixed AB state, with the additional rule that a node can- 
not change state from A to B or vice versa without going 
through the AB state. In studies of dynamics of lan- 
guage competition, the voter model gives a microscopic 
version [15] of the Abrams-Strogatz [16] model for the 
competition of two socially equivalent languages. In this 
context the third state of non-excluding options of the 
AB model is naturally associated with bilingualism [17]. 
More generally the AB model describes competition of 
two equivalent social norms which can coexist at the in- 
dividual level. 

In both models, an agent changes its state with a prob- 
ability which depends on the states of its neighbors. The 
fraction of first neighbors in state A [B, AB] of an agent is 
called the local density of A, a a [o'b-, cab]- For the voter 
model, the state AB is not allowed and the probabilities 
of a node changing state are defined as follows: 

PA^B = O'B, PB^A = CFA ■ (1) 

The AB model is defined by the following update rules: 
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In our simulations we start from random initial con- 
ditions for the state of the agents in a network with N 
nodes (see below) and we use random asynchronous node 
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update: at each time step a single node is randomly cho- 
sen and updated according to the transition probabilities 
Eq. (1) or Eqs. (2)-(3). Wc normalize time so that every 
unit of time includes N time steps. 

A question of interest is under which conditions con- 
sensus is reached (all nodes hold the same option), and 
which is the process of emergence and growth of spatial 
domains where the nodes are in the same state ( coarsen- 
ing) . Both models are symmetric by interchange of A and 
B, so that reaching consensus in either of these two states 
is a symmetry breaking process. To describe the dynam- 
ics of the system we use as order parameter the interface 
density p, which is defined as the fraction of links which 
connect nodes in different states. The ensemble average 
interface density (p) is considered, where the ensemble 
average, indicated as (•), denotes average over realiza- 
tions of the stochastic dynamics starting from different 
random initial conditions. Interface density decreases as 
domains grow in size. If one of the states becomes dom- 
inant, the interface density decreases along with the dis- 
appearing state. Zero interface density indicates that an 
absorbing state, consensus, has been reached. Coarsen- 
ing in the voter model is driven by interfacial noise, while 
for the AB model earlier results indicate that coarsening 
is curvature driven: boundaries tend to straighten out, 
reducing curvature and leading to the growth of spatial 
domains [14]. It turns out that domains of AB agents 
are never formed. Instead, AB agents place themselves 
in the interface between A and B domains. 

The dependence of the voter model dynamics on net- 
work dimensionality, disorder and degree distribution has 
been carefully studied [13, 18-20]. A main result is that 
d = 2 is the critical dimensionality for this model. This 
means that for d < 2 there is coarsening, i.e. unbounded 
growth (in the thermodynamic limit) of domains in which 
all nodes are in the same state. However for d > 2 there is 
no coarsening beyond an initial transient. In finite net- 
works of d > 2 there exist long-lived metastable states 
in which p takes a plateau value. The inverse of this 
plateau value is the characteristic size of coexisting A 
and B domains. Eventually a finite size fluctuation takes 
the system to one of the two consensus absorbing states. 
We note that complex networks are typically high di- 
mensional structures for which these metastable states 
naturally occur [13]. 

Coarsening processes leading to consensus often come 
to a halt due to the appearance of metastable states that 
can be of different nature. Coarsening and metastable 
properties depend on the dynamical model as well as on 
network characteristics. The type of mctastability en- 
countered for the voter model is characterized by the fact 
that all realizations of the process are of the same class 
(qualitatively similar) and that the metastable states 
have a finite lifetime for a finite system. For the voter 
model the mean lifetime of these states scales as t 
[13]. We call this type of metastable states dynami- 
cal metastable states. A different type of metastability, 
which we call trapped metastable states., occurs in situ- 




FIG. 1: Growth process of the network. The new vertex v 
links to one or more randomly chosen initial contacts (here 
and possibly to some of their neighbors (here k,l). 

ations in which different realizations of the process are 
of different type. While some of them follow a coarsen- 
ing process until finite size effects come into play, others 
get stuck in topological traps. The latter correspond to 
trapped metastable states that can be of two types: they 
might have a finite lifetime in finite systems, as it oc- 
curs for the AB model with stripe-like configurations in 
regular two dimensional lattices [14], or they might be 
infinitely long lived as it occurs in zero temperature ki- 
netic Ising models [21-24]. In summary, different forms 
of metastability can appear for the voter and AB models 
considered here, but every realization is expected to have 
a finite lifetime in a finite system. 



III. A CLASS OF SOCIAL TYPE NETWORKS 

Several models have been designed to capture some of 
the characteristics of social networks, based on mecha- 
nisms such as geographical proximity [8] , social similarity 
[3, 7], and local search [5, 6, 9]. A combination of random 
attachment with local search for new contacts has proved 
fruitful in generating cohesive structures as well as well- 
known features of social networks, such as assortativity, 
broad degree distributions, and community structure [9]. 
The term "community" is typically used in the context of 
groups of nodes with dense internal and sparse external 
connections; exact definitions differ [25-28]. The com- 
munity structure leads naturally to high values of the 
clustering coefficient and to positive degree-degree corre- 
lations. 

The algorithm to generate this class of networks con- 
sists of two growth processes: 1) random attachment, 
and 2) implicit preferential attachment resulting from 
following edges from the randomly chosen initial con- 
tacts. The local nature of the second process gives rise to 
high clustering, assortativity and community structure. 
Starting from any small connected seed network of iVo 
vertices, new nodes are added as follows (see Fig. 1): i) 
Pick Uinit > 1 random nodes as initial contacts, ii) Pick 
Usee > neighbors of each initial contact as secondary 
contacts, iii) Connect the new node to the initial and 
secondary contacts. 

Throughout this paper, we will use the standard pa- 
rameters [9]: the number of initial contacts is selected 
according to the probabilities p{ninit = 1) = 0.95, 
p{ninit — 2.) — 0.05; and the number of secondary con- 
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FIG. 2: Left: A partial view of the network centered on a 
randomized selected node. Right: Average number {n{s)} of 
fc-clique-communities of size s for fc = 4 (□) and k — 5 (0)i in 
networks of size N = 10 000, averaged over 400 realizations. 

tacts from each initial contact, risec, is chosen from a 
uniform probabihty distribution between and 3; the 
initial seed contains A^o = 10 nodes. 

The degree distributions of the resulting networks are 
found to decay slower than exponential [9] . Using the k- 
clique algorithm [28] for detecting communities, a broad 
distribution of community sizes is found in the model 
(Fig. 2). 

For reference, we use randomized versions of the same 
networks, where the degree sequence is kept intact but 
edges are randomly rewired under the restriction that 
the network must stay connected [29]. This eliminates 
community structure, clustering, and degree correlations. 
The randomized networks are therefore locally treelike 
with very few loops. 



IV. RESULTS 

We have considered the update rules Eqs. (1) for the 
voter model, or Eqs. (2)- (3) for the AB model in the 
class of networks described above. We followed the de- 
velopment over time of the interface density and of the 
fraction of runs that had not yet reached consensus at 
any particular time. When results for the original and 
randomized networks differ, we can conclude that struc- 
tural characteristics other than the degree distribution 
are responsible for the differences. 




FIG. 3: Time evolution of the average interface density in 
networks with communities (solid symbols) and randomized 
networks (empty symbols) with the same degree sequences, 
(a) Voter model. Network sizes increase from left to right: 
iV = 100, 400, 2500, 10000. Averages are taken over 100 
different realizations of the model network, with 10 runs in 
each. Inset: time to reach consensus scales with network size 
as r ~ A'^'^, 7 « 0.96 for the randomized and 7 « 0.98 for the 
original networks, (b) AB model. Network sizes increase from 
left to right: iV = 100, 400, 2500, 10000, 40000. Averages 
taken over 400 — 5000 realizations (depending on system size) 
of the model network, and with 10 runs in each. 



A. Interface density 

The average interface density (p) on the class of net- 
works considered here, and on their randomized counter- 
parts is shown in Fig. 3. For the voter model (Fig. 3a), 
we obtain that the structure of the network does not 
alter the qualitative behavior. In both classes of net- 
works we observe plateau values of (p) associated with 
dynamical metastable states. Still, the plateau value for 
networks with community structure is lower than for the 
randomized networks, indicating that the typical size of 
spatial domains where agents are in the same state is 



larger. We also observe in both cases that finite size fluc- 
tuations drive the system to an absorbing state. The 
characteristic time to reach consensus (mean lifetime of 
the metastable state) depends on network size but it does 
not depend sensitively on network structure. The inset in 
Fig. 3a shows that the time to reach consensus depends 
linearly on network size for networks with communities 
and their randomized counterparts [30]. These results 
support the earlier finding made on networks without 
mesoscopic structure that effective dimensionality domi- 
nates voter model behavior [13]. 

Figure 3b shows the average interface density for the 
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AB dynamics. We observe significant differences be- 
tween tlie original and the randomized version networks: 
a plateau value of (p) is observed for randomized net- 
works, while a first dynamical stage of coarsening where 
spatial domains grow in size is found for large networks 
with communities. The plateau observed in randomized 
networks indicates that a dynamical metastable state of 
the class found in the voter model for both types of net- 
works is rapidly reached. Moreover, in the randomized 
networks there is a fast decay towards an absorbing state 
with a characteristic time to reach consensus almost in- 
dependent of system size. For the networks with a com- 
munity structure we observe two dynamical stages in the 
evolution of (p). After an initial power law associated 
with coarsening there appears a second power law tail in 
the approach to the absorbing state. This last power law 
decay indicates that the mean lifetime to reach consen- 
sus for the AB model docs not characterize the dynamics 
on these networks and that metastable states exist at all 
time scales, as we discuss below. Additionally, the dif- 
ference with the randomized networks in several orders 
of magnitude for the extinction times, which increases 
with system size, shows that the network with communi- 
ties slows down the dynamics significantly. All together 
these results manifest a sensitivity of the AB dynamics 
to the mesoscopic network structure which is not found 
for the voter dynamics. 



B. Fraction of alive runs 
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FIG. 4: Fraction of alive runs in time for networks with com- 
munities (solid symbols) and randomized networks (empty 
symbols). AB model (double logarithmic plot); system sizes 
iV = 100, 400, 2500, 10000 from left to right, with averages 
taken over different realizations of the network (400 — 5000 
depending on system size), with 10 runs in each. Inset: voter 
model (semilogarithmic plot). System sizes A*' = 2500, 10000. 
Averages are taken over 100 different realizations of the net- 
works, with 10 runs in each. 



level. Such structure seems to give rise to a number of 
traps that cause trapped metastable states at all time 
scales. To substantiate this claim we next look at some 
detailed dynamics. 



Figure 4 shows the fraction P{t) of realizations still 
alive at time t, i.e. the fraction of realizations which 
have not reached the absorbing state. For the voter 
model, the fraction of alive runs decreases exponentially 
in both the original and randomized networks (Fig. 4- 
inset), in agreement with previous results for the voter 
model in high dimensional complex networks [13]. A 
rather different result is found for the AB model (Fig. 4). 
In our class of networks, we find a power law behavior 
P{t) ~ t^", a « 1.3, so that a mean lifetime of the 
realizations of the AB dynamics does not give a charac- 
teristic time scale. At any time there are live realizations 
which have not reached the absorbing state. Different 
parametrizations of the network model (not shown) pro- 
duce the same qualitative phenomenon: we have modi- 
fied the number of secondary contacts from each initial 
contact, risecj using uniform probability distributions be- 
tween and 1,2,4, obtaining also a power law of the 
distribution of alive runs with an exponent smaller than 
2, which indicates the robustness of this result. This 
behavior is different from the usual exponential decay of 
the tails of P{t) observed for the voter, and AB dynamics 
either in regular, small world [14], random or Barabasi- 
Albert scale-free networks (not shown), and reflects the 
existence of metastable states at all time scales. This 
fact indicates that the anomalous lifetime distribution is 
linked to the structure of the network at a mesoscopic 



C. Discussion 

Further understanding of the dynamical process can 
be obtained by considering the measure called overlap, 
O [4]. This characteristic of a link between two nodes 
tells us essentially which fraction of their neighbors is 
shared by the nodes. Within a community, nodes tend 
to share many neighbors, and thus overlap is high, while 
edges between communities will have low or zero over- 
lap. Considering dynamics of competing options on a 
network, the overlap can be used to identify spatially ho- 
mogenous domains in the network: if the average overlap 
(O) of the links in the interface between domains is low, 
we may assume that the domain boundaries follow the 
community boundaries. On the other hand, if the over- 
lap at the interfaces is high, it indicates that nodes within 
communities are in different states. For the voter model 
dynamics we have found that the average overlap of inter- 
face links drops to about 80 percent of the average value 
(O) = 0.27 of the whole network, while in the AB model 
it drops to under 70 percent. This indicates that in both 
models the interfaces between domains lie preferably in 
low overlap links, so that domains of the same option fol- 
low the community structure, but in the AB model these 
domains arc correlated with the communities closer. 

The difference between the two dynamics is better un- 




FIG. 5: Snapshots of the dynamics, with nodes in state A in 
black, B in grey, and AB in white circled in black. Simula- 
tions start from random initial conditions. Left: voter model. 
Right: AB model. 



derstood by looking at snapshots of the dynamics (Fig. 5) 
which show the characteristic behavior for each of the 
models, starting from random initial conditions = 0). 
In the voter model (left) the homogeneous domains of 
nodes with the same option appear to follow the commu- 
nity structure, but a particular community (topological 
region) may change the option adopted by the commu- 
nity rather quickly (t = 50, 60, 70). At variance with this 
behavior, in the AB model (right) spatial domains grow 
and homogenize steadily in a community without much 
fluctuation. For this dynamics, communities that have 
adopted a given option, and which are poorly linked to 
the rest of the network, take a long time to be invaded by 
a different option, acting therefore as topological traps. 
As an example of this we show two long lived trapped 
metastable state at t ~ 430 and t = 1000, where the 
interface stayed relatively stable for a prolonged period 
(^ 100 and ^ 1000 time steps, respectively). These dif- 
ferent behaviors reflect in the community structure two 
different interfacial dynamics: interfacial noise driven dy- 
namics for the voter model, and curvature driven dynam- 
ics for the AB model with agents in the AB state at the 
interfaces. 




FIG. 6: (a) Time evolution for the AB model of the aver- 
age interface density on different realizations of the network 
with 2500 agents; 20000 runs on each (empty symbols). The 
extreme cases were selected as examples of networks where 
trapped metastable states (see text) are found often (0)i s^nd 
found rarely (A). For comparison, the average over 500 net- 
works (10 runs on each) is also shown (■). Inset: time evo- 
lution for the voter model of the average interface density for 
four realizations of the networks of 2500 agents; 5000 runs on 
each network, (b) Time evolution of the interface density in 
single realizations of the AB dynamics on a network with 2500 
agents. A class of realizations decay to the absorbing state 
after a coarsening stage (solid black lines), while others fall in 
long lived trapped metastable states. The latter display sev- 
eral plateaus, indicating hierarchical levels of ordering before 
reaching the absorbing state, or cascading between several 
trapped metastable states. 



Different realizations of the algorithm to construct the 
social type network produce different detailed structures 
of the network. The power-law for the fraction of alive 
runs in Fig. 4 is a statistical effect of the average over such 
realizations. The time evolution of the average interface 
density on single realizations of the network, (p) , is shown 
for the AB dynamics in Fig. 6a. We observe different be- 
haviors in the second stage of the decay of (p) depending 
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on the specific realization of the network: from broad 
tails to exponential- like decays, with an intermediate be- 
havior. On the other hand, and in agreement with our 
previous discussion, the voter model dynamics (Fig. 6a, 
inset) is not sensitive to the details of the network struc- 
ture. For the AB model some realizations of the net- 
work topology produce particularly long lived mctastable 
states, while in others, corresponding to exponential-like 
decay of (p) , they are observed rarely. Plots of the inter- 
face density of individual runs on a given network show 
a class of realizations with different plateaus (ordering 
levels) where the system gets trapped for a long time 
(Fig. 6b). These trapped metastable states, analogous to 
those displayed in Fig. 5-right, correspond to the struc- 
ture in the network. The variety of traps and associated 
different lifetimes seems to be the mechanism that causes 
an anomalous power law distribution for the lifetimes. 

We note that although the details of each network real- 
ization matter for the occurrence of trapped metastable 
states, the community size distribution detected by the 
fc-clique-pcrcolation method [28] is the same for all the 
network realizations that we have considered. This and 
other available statistical methods seem not to be suffi- 
cient to discern between the network topologies produc- 
ing many or few trapped metastable states. 

V. SUMMARY AND CONCLUSIONS 

We have considered two dynamical models, the voter 
and the AB model, in order to study metastable states 
and the role of community structure in the dynamics of 
consensus processes. The voter model dynamics, driven 
by interfacial noise, is not particularly sensitive to the 
mesoscopic structure of the network: we find that all re- 
alizations of the dynamics are of the same class, leading 
to a type of dynamical metastable states shared by other 



complex networks of high dimensionality without degree 
correlations. On the contrary, for the AB dynamics we 
find different classes of realizations leading to a power 
law distribution for the times to reach consensus. This 
is explained in terms of trapped metastable states as- 
sociated with the structure of the network. Our result 
implies that a mean lifetime for these states docs not 
give a characteristic time scale of the ordering dynam- 
ics. We note that a mean lifetime does not exist for the 
zero-temperature kinetic Ising model dynamics on reg- 
ular or complex networks [23], due to realizations that 
lead to trapped metastable states of infinite lifetime in 
finite systems. The novelty of our finding is that we have 
realizations with any lifetime. For the AB model in a 
regular 2D lattice trapped mctastable states with stripe- 
like configuration have been found [14], but in that case 
the distribution of lifetimes is exponential: P{t) ^ e~"* 
and the mean lifetime is representative of the dynamics. 
The power-law distribution for the lifetimes originates 
here in the multiplicity of different traps that reflects the 
mesoscopic structure of the networks. Simpler configu- 
rations of community structure should be considered in 
the future in order to gain a deeper understanding of the 
microscopic mechanisms underlying consensus dynamics. 
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